A Cross-correlation method to search for 
gravitational wave bursts with AURIGA and Virgo 

THE AURIGA COLLABORATION 

M. Bignotto 1 2 , M. Bonaldi 34 , M. Camarda 5 , M. Cerdonio 1 2 , 
L. Conti 12 , M. Drago 12 , P. Falferi 3 - 4 , N. Liguori 12 , S. Longo 6 , 
R. Mezzena 47 , A. Mion 4 7 , A. Ortolan 6 , G.A. Prodi 4 7 , V. Re 4 7 , 
F. Salemi 4 * 7 , L. Taffarello 2 , G. Vedovato 2 , A. Vinante 3 - 4 , S. 
Vitale 47 , J.P. Zendri 2 

1 Dipartimento di Fisica, Universita di Padova, Via Marzolo 8, 35131 Padova, Italy 

2 INFN, Sezione di Padova, Via Marzolo 8, 35131 Padova, Italy 

3 Istituto di Fotonica e Nanotecnologie, CNR-Fondazione Bruno Kessler, 1-38050 
Povo (Trcnto), Italy 

4 INFN, Gruppo Collegato di Trcnto, Sezione di Padova, 1-38050 Povo, Trcnto, Italy 

5 Dipartimento di Ingegneria Informatica, Universita di Padova, Via G. Gradenigo 
6a, 35131 Padova, Italy 

6 INFN, Laboratori Nazionali di Legnaro, 35020 Legnaro, Padova, Italy 

7 Dipartimento di Fisica, Universita di Trento, 1-38050 Povo, Trcnto, Italy 

THE Virgo COLLABORATION 

F.Acernese 7 ' 9 , M.Alshourbagy 1516 , P.Amico 1314 , F.Antonucci 19 , 
S.Aoudia 10 , P.Astone 19 , S.Avino 7 ' 8 , L.Baggio 1 , G.Ballardin 2 , 
F.Barone 79 , L.Barsotti 15 16 , M.Barsuglia 11 , Th.S. Bauer 21 , 
S.Bigotta 15 16 , S.Birindelli 1516 , C.Boccara 12 , F.Bondu 10 , L.Bosi 13 , 
S.Braccini 15 , C.Bradaschia 15 , A.Brillet 10 , V.Brisson 11 , 

D. Buskulic 1 , G.Cagnoli 3 , E.Calloni 78 , E.Campagna 3 ' 5 , 
F.Carbognani 2 , F.Cavalier 11 , R.Cavalieri 2 , G.Cella 15 , 

E. Cesarini 34 , E.Chassande-Mottin 10 , A.-C.Clapson 11 , F.Cleva 10 , 
E.Coccia 23 ' 24 , C.Corda 1516 , A.Corsi 19 , F.Cottone 1314 , 
J.-PCoulon 10 , E.Cuoco 2 , S.D'Antonio 23 , A.Dari 13 14 , V.Dattilo 2 , 
M.Davier 11 , R.De Rosa 78 , M.Del Prete 15 17 , L.Di Fiore 7 , A.Di 
Lieto 1516 , M.Di Paolo Emilio 23 ' 25 , A.Di Virgilio 15 , M.Evans 2 , 
V.Fafone 23 24 , I.Ferrante 15 16 , F.Fidecaro 15 16 , I.Fiori 2 , 
R.Flaminio 6 , J.-D.Fournier 10 , S.Frasca 19,2 °, F.Frasconi 15 , 
L.Gammaitoni 1314 , F.Garufi 7 ' 8 , E.Genin 2 , A.Gennai 15 , 
A.Giazotto 2 15 , L.Giordano 78 , V.Granata 1 , C.Greverie 10 , 
D.Grosjean 1 , G.Guidi 3,5 , S.Hamdani 2 , S.Hebri 2 , H.Heitmann 10 , 
P.Hello 11 , D.Huet 2 , S.Kreckelbergh 11 , P.La Penna 2 , M.Laval 10 , 



A Cross- correlation method to search for gravitational wave bursts with AURIGA and Virgo2 



N.Leroy n , N.Letendre 1 , B.Lopez 2 , M.Lorenzini 3 ' 4 , V.Loriette 12 , 
G.Losurdo 3 , J.-M.Mackowski 6 , E.Major ana 19 , C.N.Man 10 , 
M.Mantovani 1716 , F.Marchesoni 1314 , F.Marion 1 , J. Marque 2 , 
F.Martelli 3 ' 5 , A.Masserot 1 , F.Menzinger 2 , L.Milano 7 ' 8 , 
Y.Minenkov 23 , C.Moins 2 , J.Moreau 12 , N.Morgado 6 , S.Mosca 78 , 
B.Mours 1 , I.Neri 1314 , F.Nocera 2 , G.Pagliaroli 23 , C.Palomba 19 , 

F. Paoletti 2 15 , S.Pardi 78 , A.Pasqualetti 2 , R.Passaquieti 1516 , 
D.Passuello 15 , F.Piergiovanni 3 ' 5 , L.Pinard 6 , R.Poggiani 1516 , 
M.Punturo 13 , P.Puppo 19 , P.Rapagnani 19 20 , T.Regimbau 10 , 
A.Remillieux 6 , F.Ricci 19,2 °, I.Ricciardi 7,8 , A.Rocchi 23 , 
L.Rolland 1 , R.Romano 7 ' 9 , P.Ruggi 2 , G.Russo 7 ' 8 , S.Solimeno 7 ' 8 , 
A.Spallicci 10 , B.L.Swinkels 2 , M.Tarallo 1516 , R.Terenzi 23 , 
A.Toncelli 1516 , M.Tonelli 1516 , E.Tournefier 1 , F.Travasso 1314 , 

G. Vajente 1816 , J.F.J, van den Brand 21 ' 22 , S. van der Putten 21 , 
D.Verkindt 1 , F.Vetrano 35 , A.Vicere 35 , J.-Y.Vinet 10 , H.Vocca 13 , 
M.Yvert 1 

^aboratoire d'Annccy-lc-Vicux dc Physique des Particules (LAPP), IN2P3/CNRS, 

Universite de Savoic, Annccy-le-Vieux, France 

2 European Gravitational Observatory (EGO), Cascina (Pi), Italia. 

3 INFN, Sezione di Firenze, Sesto Fiorentino, Italia. 

4 Universita degli Studi di Firenze, Firenze, Italia. 

5 Universita degli Studi di Urbino " Carlo Bo" , Urbino, Italia. 
6 LMA, Villeurbannc, Lyon, France. 

7 INFN, sezione di Napoli, Italia. 

8 Universita di Napoli "Federico II" Complesso Universitario di Monte S.Angclo, 
Italia. 

9 Universita di Salerno, Fisciano (Sa), Italia. 

10 Departement Artemis - Observatoire de la Cote d'Azur, BP 4229 06304 Nice, 
Cedex 4, France. 

11 LAL, Univ Paris-Sud, IN2P3/CNRS, Orsay, France. 

12 ESPCI, Paris, France. 

13 INFN, Sezione di Perugia, Italia. 

14 Universita di Perugia, Perugia, Italia. 

15 INFN, Sezione di Pisa, Italia. 

16 Universita di Pisa, Pisa, Italia. 

17 Universita di Siena, Siena, Italia. 

18 Scuola Normale Superiore, Pisa, Italia. 
19 INFN, Sezione di Roma, Italia. 
20 Universita "La Sapicnza", Roma, Italia 

21 National institute for subatomic physics, NL-1009 DB Amsterdam, The 
Netherlands. 

22 Vrije Universiteit, NL-1081 HV Amsterdam, The Netherlands. 
23 INFN, Sezione di Roma Tor Vergata, Roma, Italia. 

24 Universita di Roma Tor Vergata, Roma, Italia. 

25 Universita dellAquila, L Aquila, Italia. 



e-mail virginia.re@lnl.infn.it 



A Cross-correlation method to search for gravitational wave bursts with AURIGA and Virgo3 
Abstract. 

We present a method to search for transient GWs using a network of detectors 
with different spectral and directional sensitivities: the interferometer Virgo and the 
bar detector AURIGA. The data analysis method is based on the measurements of 
the correlated energy in the network by means of a weighted cross-correlation. To 
limit the computational load, this coherent analysis step is performed around time- 
frequency coincident triggers selected by an excess power event trigger generator 
tuned at low thresholds. The final selection of GW candidates is performed by a 
combined cut on the correlated energy and on the significance as measured by the 
event trigger generator. The method has been tested on one day of data of AURIGA 
and Virgo during September 2005. The outcomes are compared to the results of a 
stand-alone time-frequency coincidence search. We discuss the advantages and the 
limits of this approach, in view of a possible future joint search between AURIGA and 
one interfcrometric detector. 



1. Introduction 

We present a study on the performances of a gravitational wave (GW) observatory 
composed by a hybrid network of detectors. In particular, we focus on the possibility to 
use a resonant detector to perform GW observations with one interferometric detector. 
This could be an opportunity in the scenario after LIGO S5 run and the first Virgo 
science run, when most of the interferometers will be shut down for upgrading: current 
plans are that GEO will be kept in operation till the start of the LIGO S6 and the 
second Virgo science runs, supported by LIGO Hanford 2k detector over weekends. In 
this sense, we present a case study on joint observations between AURIGA and Virgo 
on a test period of 24 hrs. 

In the past years, various searches for GW signals have been independently 
performed by networks of resonant bars [TJ [21 E] or interferometers jU E] . There have been 
also some attempts to perform burst searches among detectors with different spectral 
sensitivity and orientation: by TAMA and the LIGO Scientific Collaboration (LSC) [S], 
by AURIGA and the LSC [3 El [9] and by the INFN bars and the Virgo Collaboration 

mi- 

The proposed network search strategy takes as a starting point the Wave- 
Burst +CorrPower [121 H3] search used by LSC for the S3 and S4 analyses [HI E]- That 
search was greatly innovative: a two-step search composed of an ExcessPower-like event 
trigger generator plus a cross-correlation test which allowed an efficient reduction of 
false alarms. In that case, however, the detectors partecipating to the network were 
almost aligned and had a similar spectral sensitivity. An extension of such methodology 
to the case of similar but misaligned detectors has been discussed in literature [T5] . 
The novelty of our work consists in a further generalization to detectors with different 
spectral sensitivities, so that it can be implemented between a resonant bar and an 
interferometer. To better characterize the method, we compare its performances with 
those of a simple time-frequency coincidence search. 
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The paper is organized as follows: in section 2 we introduce the search method. 
Section 3 presents an overview of the exchanged data and summarizes the main steps 
of the network pipeline and of the tuning performed on chosen test-statistics. Results 
and conclusions are presented in section 4 and 5 respectively. 

2. The search method 

The G W search method described in this paper is characterized by two main parts: 
the event trigger generator, whose role is to select a subset of interesting triggers and 
a coherent analysis. The trigger search is based on Waveburst [16], an excess power 
algorithm based on the wavelet decomposition in the time-frequency plane. In the 
present work, Waveburst has been used in coincidence mode, i.e. the algorithm selects 
time-frequency coincident excesses of power between the two detectors. The step of 
coherent analysis is based on a cross-correlation test between data streams weighted by 
a combination of the strain sensitivities of the two detectors (XCorr). 

Our method assumes that the GW components at earth can be parametrized as 

h+(t) = ho(t) ■ ca#(t)] K (t) = e ■ h (t) ■ sm[^(t)} (1) 

where h (t) and ip(t) are the time- varying amplitude and time- varying phase, common 
to both polarization components, and e is the ratio of the cross and plus amplitudes. 
A large class of GW signals can be parametrized as described above, including linearly, 
elliptically and circularly polarized GWs, even with sweeping frequencies. 

The strain produced on the detector a by an incoming burst signal with polarization 
components in the wavefront frame h +<x (t) is: 

h a (t) = F a+ ■ h + (t) + F ax ■ h x (t) (2) 

where F a+ and F ax are the antenna pattern functions [10] (and references therein). 
Following [T7], eq|2] becomes: 

K{t) = h {t) ■ R a (6, 0, e) • cos^(t) + f Q (0, 0, e)] (3) 

where R a is a directional sensitivity, £ Q is a phase shift and (6, 0) is the location of 
the GW source [|j The reconstructed strain at the input of two detectors, a and j3, 
is x a ,p{t) = ha,p(t) + n a fi(t), where n a ,/3 are the two independent noises. It has been 
shown in [TH] that the following linear combination of the two reconstructed strains, 
called null stream, cancels the signal: x nu u(t) = x a {t)Rp — xp(t + f)R a , where t' 
includes the light travel time and a suitable fraction of the typical GW period, so 
that + C« = + t') + 6/i 

| With some algebra, it is easy to find that R a = (F a+ ) 2 + (F ax ■ e) 2 and £ Q = —arctan[ F p X + e ]. 
§ An additional assumption is required here: the envelope of the GW signal be smooth in time so 
that h (t) ~ h (t + t'). Whenever one considers a cross-correlation with a narrow-band detector, this 
approximation is automatically verified since the reconstructed strain at input has to be bandlimited 
by a suitable filter (see fig(T|. 
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We introduce the Fourier transform of the null stream and normalize it to its 
variance per unit bandwidth, a 2 ull (u): we properly filter x a (t) obtaining the weighted 
reconstructed strain at detector a: 

x aw {u) = AM = W (4 ) 

where S aj p are the noise power spectral densities of the detectors in terms of GW strain. 
Hence, the normalized null stream is x nu u )W (t) = x atW (t)Rp — xp >w (t + t')R a . 

One well known method to search for GW signals in the data relies on the 
minimization of the null energy [T7J HE], E n uii = J dtx\ ull w (t) , where the time integral 
is performed on the signal duration plus any typical response time of the narrower band 
detector. The null energy can be expressed in terms of the correlated energy and the 
incoherent energy of the network: E nu u = —E cor +Ei nc . The former is the contribution of 
the cross-correlation of the detectors, E cor = 2 / dtx ayW (t)xp jW (t + t')R a Rp. The latter is 
the auto-correlation contribution of the detectors, E inc = J dt(x 2 a (t)R 2 p+x 2 Jt+t')R 2 a ). As 
discussed in [18], a GW candidate is selected against background events more efficiently 
by maximizing E cor rather than by minimizing E nu u. In fact, E nu u can take low values 
even for accidental events with small E cor and E inc \ instead, for detectable GW signals, 
we expect a higher E inc , almost balanced by a positive E cor . For these reasons, this 
coherent step of network analysis is based on the maximization of the correlated energy 
E cor in our null stream. 

In principle, E cor depends on 9, (ft and e of the source through t', R a and Rp. 
However, we checked that in the case of random polarized GW radiation, emitted 
by sources distributed either in the galaxy or uniformly in the sky, we can follow an 
approximated maximization procedure of E cor assuming R a ~ Rp. 

The main reason is that AURIGA is limiting the common bandwidth and at the 
same time shows a worse strain power spectral density; therefore, the a nu u{uj) defined in 
eqjijis dominated by the strain noise of AURIGA unless Rau/Rv 3> 1- This condition 
however is occurring rarely enough: Rau/Rv < 2.4 with 90 % probability for a galactic 
population of sources, and very similar values occur also considering the LIGOs or 
GEO interferometers in place of Virgo. Figjl] shows the a nu u(uj), corresponding to 
Rau/Rv — {0, 1,2.4} . These curves are quite similar and therefore it is advisable to 
approximate R a = Rp in eq|4j i.e. to consider filtered data streams: 

x'^wiu) = , (5) 
'S a (uj) + Sp(u) 



to compute the cross-correlation. In addition, to be robust against noise fluctuations and 
relative calibration errors between AURIGA and Virgo, we considered the r-statistics of 
such filtered data streams: 

fdtx> aiW (t)x^ w (t + t>) 

r a P = / = (O) 

dTX* w {t) J dTX% w {t) 

where we integrate over a time window of fixed duration and search for local maxima 
as a function of the time t and time shift t! between the data streams. 
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Figure 1 : Strain noise power spectral densities of AURIGA and Virgo detectors on September 
the 15 , 2005 (black continuous lines). The gray curve shows the filter <J nu u(uj) with the 
approximation Rau = Rv- The exact shape of the spectral filter is lower limited by the 
AURIGA strain noise, corresponding to sources with Rau Rv- The light gray curve shows 
the exact spectral filter for sources with Rau = 2.4i?y and can be interpreted as an upper 
limiting curve since the probability that Rau < 2ARy is 90% for a source distributed in the 
Milky Way with random polarization. We conclude that the spectral filter with Rau = Rv is 
a satisfactory approximation for our purposes. 

3. Network Analysis 

We analysed the 24 hours of data starting from the UTC time 14 Sep 2005 23:11:27. 
In this time period, Virgo was performing its seventh commissioning run (C7) while 
AURIGA was in stable operation. This affects mainly the quality of the data of Virgo 
for the presence of periods of instrumental unlock reducing the observation time from 
24 to 16 hours and 33 mins. Moreover, the Virgo Collaboration provided a list of vetoes 
to flag the most noisy periods, corresponding to 14.5% reduction of the live time. The 
Shh sensitivity curves for the two detectors in the [850 -j- 950] Hz band on September 
lh th are shown in figfl] 

3.1. Pipeline description 

Fig|2] shows the scheme of the implemented pipeline: as a first step, the raw data from 
the two detectors are whitened and bandpassed to match the AURIGA bandwidth (i.e. 
850 - 950 Hz), then they are fed in as input to WaveBurst [Tfi"] . 

WaveBurst has been configured to search for coincident excesses of power on the 
two data sets over 4 wavelet decomposition levels, corresponding to time /frequency 
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resolutions ranging from 4.916 ms by 101.72 Hz to 39.32 ms by 12.72 Hz; it then produces 
the list of triggers, containing all the trigger parameters (e.g. the central time in GPS 
seconds, the central frequency, the frequency range, the geometric significance, the 
Signal to Noise Ratio, SNR, etc. ); for each of those events, the last step of the pipeline, 
XCorr, selects the corresponding chuncks of weighted filtered data (see sec. [2J and, by 
sliding on t and t' the two time series [JJ it searches the maximum of the cross-correlation 
coefficient, r, computed over a given time window. 

Within our approach, Z and r are the main test statistics which characterize each 
event at the output of the pipeline. A combined threshold on these quantities sets the 
overall false alarm rate (FAR) and detection efficiency of the network. Other internal 
parameters have been set in order to have neglegible impact on the final FAR and 
detection efficiency. 
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Whitened Date 



Filtered Data 
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Event Trigger 
Generator 
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Figure 2: Block diagram of the pipeline. The two detectors whitened data sets are passed to 
WaveBurst that produces a list of coincident triggers, each characterized by the geometrical 
significance, Z. The XCorr takes the two weighted filtered data chuncks around the coincident 
trigger times and slides them looking for the maximum of the cross-correlation coefficient, r. 



3.2. Background and detection efficiency estimation 

The background is estimated by applying the pipeline to 4000 time-shifted data sets 
with a total live-time of 2760 days, with shifts ranging from -2123s to +2 123s in steps 
of m 1 s. This choice allows to have a sufficiently large number of resamplings, while 
preserving the main characteristics of the data sets. 

The detection efficiency has been computed over a number of different waveforms 
(mainly elliptically polarized damped sinusoids) by means of software injections: the 
signal central frequency /o ranges within the bar bandwidth (850-950 Hz) and the decay 
time r spans at most a few tens of milliseconds, with random inclination and polarization 

|| For the definition of the significance, see eq.s 10 and 11 in [T5] 

% Between Virgo and AURIGA, the light travel time is < 0.8 ms and the phase shift is < 0.5 ms; to 
be conservative we set the relative time slide to \t'\ < 1.5 ms. 



A Cross-correlation method to search for gravitational wave bursts with AURIGA and Virgo8 

angles. We have considered source populations distributed at the Galactic Center, in 
our galaxy according to its mass distribution [19] and uniform distributed in the sky. 

For each waveform, we have created two sets of 994 signals, approximately spaced 
by 60 s: the tuning and the estimation sets. 

3. 3. Tuning of the analysis 

The aim of the tuning phase is to optimize the pipeline parameters in order to achieve a 
low FAR, while preserving the detection efficiency The search pipeline has been tuned 
using 2000 time-shifts randomly selected out of the initial 4000 and the tuning set of 
injections. The parameters set in this tuning phase have then been used on the second 
half of the time shifts to get the final FAR and on the estimation set of injections to 
calculate the network efficiency 

Cross- correlation Window Among the pipeline parameters, the cross-correlation 
window is one of the most important. The general view is that a short time window 
may be unsuitable for a long signal, as a significant part of its power may be cut off by 
the window itself. At the same time, if the integration time span is too large, the signal 
may be diluted in it. Indeed, the weighted cross-correlation filter induces a correlation 
between the 2 data series over a time scale of the order ~ 30ms. Therefore, even for 
delta-like signals, we need a cross-correlation window of at least 60 ms. We have set it 
to 100 ms after testing the performances of various time windows ([25, 50, 100, 200, 400, 
800]ms) in terms of FAR and efficiency for narrow and wide band simulated signals. 

Combined threshold We set the FAR to ~ 2/yr: given our short observation time, we 
could not target a smaller rate, since only 7 background events were surviving. We have 
not applied the Virgo list of vetoes, as from tests performed on the tuning set of the 
shifts, we measured a ~ 30% FAR reduction only, while losing 14.5% of the live time. 

The target FAR has been achieved by setting a combined (Z, r) threshold on the 
event list selected by the pipeline. In fig|3j we show an example of how a population 
of injected signals (squares) detected by the pipeline is distributed with respect to the 
background events (dots) for two waveforms. The dashed line shows the optimized 
combined threshold on (Z,r): such threshold is defined by a simple linear equation: 
Z = mr + q, where the parameters m, q are set in such a way as to achieve the maximum 
efficiency at our target FAR. We underline here that, since AURIGA performance over 
the short time period of the selected data is quite stable, its noise distribution does not 
change significantly over different time spans and hence the choice of our threshold is 
fairly general. 

4. Results 

By applying the tuned analysis over the estimation set of the time-shifted resamples, 
we get a final FAR = 1.6/yr with a statistical sigma of 0.6/yr. The detection efficiency 



A Cross- correlation method to search for gravitational wave bursts with AURIGA and Virgo9 




Figure 3: Geometrical significance versus r-statistic for background events (dots) and two 
injected Damped Sinusoids (squares), a): DS / = 914 Hz, r = 1 ms, h rss = 1 x 10~ 19 1/y/Hz. 
b): DS / = 930 Hz, r = 30 ms, h rss = 1 x 10~ 19 1/VHz. Circled dots are background 
events resulting above threshold (dashed line) and hence contributing to our FAR. Asterisks 
are instead injected events not passing the threshold. 



of our network is shown in figjl^i, for the 3 injected source populations of damped 
sinusoids waveforms. The 50% network detection efficiency is achieved in the range 
h™f e [ 3 x lCT 20 , 1.3 x 10~ 19 ]l/VHz- The efficiency is slightly dependent on the 
source populations, while it strongly depends on the signal duration, as expected due 
to the AURIGA narrower bandwidth. 

4-1. Comparison with the time-frequency coincidence pipeline 

We have compared the results of this analysis (i.e. detection efficiency, FAR, observation 
time) with those of a simple time-frequency coincidence of single detector triggers. In the 
latter pipeline, we have analysed the two data streams separately using WaveBurst as 
event trigger generator and the resulting triggers retained if coincident in time (At = 64 
ms) and frequency (A/ = 25Hz). We have then applied the list of vetoes for Virgo C7 
data to the surviving triggers, reducing by a factor « 5 the accidental coincidences at 
the cost of 14.5% of live time. Finally, an optimized threshold on SNR has been applied 
on the resulting events from the tuning set, in order to get the target FAR~ 2/yr. The 
best efficiency has been achieved by setting the threshold at SNR > 5.6 for AURIGA 
and SNR > 14.5 for Virgo. On the estimation sets, the FAR results 5.3 ± 1.3/yr 
and the corresponding detection efficiencies are reported on in figj4jb for a Galactic 
Center source population. At similar FAR, the proposed pipeline outperforms the time- 
frequency coincidence for the tested waveforms in detection efficiency 
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Figure 4: a): Detection efficiencies of the WaveBurst+XCorr pipeline vs h rss for Galactic 
Center, all sky and Galactic Disk source populations and for the two signals (930 Hz, 30ms 
and 914 Hz, 1ms). On the upper x-axis, the h rss is converted to energy emitted in GWs at 
Galactic Center, using eq.B6 in [IT]. The network response is comparable for the 3 populations, 
but it is strongly dependent on the signal duration, b): Comparison between the detection 
efficiencies for the 2 pipelines at FAR ~ 2/yr. Solid line: WaveBurst+XCorr pipeline with 
a combined threshold on (Z,r). Dashed line: time-frequency coincidence search with two 
separate thresholds on the SNR of the detectors. For the time- frequency search only, the 
Virgo vetoes were necessary to reduce the FAR, while the WaveBurst+XCorr search results to 
be robust enough to reduce the FAR also in presence of the extra noise associated to the vetoed 
time periods. The plotted efficiencies have been normalized to the corresponding live-times, 
accounting for the loss due to the vetoes. 



5. Conclusions 

We investigated the performances of a joint search for GW bursts by a narrow-band 
resonant detector and a wide-band interferometer using one day of data taken by 
AURIGA and Virgo in September 2005. The data analysis method based on cross- 
correlation outperforms a simpler time-frequency coincidence search: it achieves better 
efficiencies at equal false alarm rate, mostly because it is by far more powerful in 
discriminating accidental coincidences. 

This search monitors at 50% efficiency galactic sources emitting ^7x 10~ 4 +■ 7 x 
lO _3 M c 2 in GW bursts of 30 +■ 1 ms decay time, provided that their strongest Fourier 
components are in the AURIGA bandwidth. The resulting detection efficiency of the 
hybrid network is limited by the less sensitive detector, AURIGA, to 5 + 10 times larger 
amplitudes with respect to a Virgo only search on a larger C7 data set [2T] for pulses 
of ~ 1ms duration. On the other hand, there is a clear advantage of this hybrid search 
in the background reduction with respect to a single detector search, as it allows to 
identify candidate events with high statistical confidence. 

The main limitations of this methodological study comes from the short duration 
of the data set used, which prevented us from investigating false alarm rates lower than 
a few per year. We expect that the efficiency of this methodology would take only a 
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small benefit from the achieved progresses in sensitivity of the interferometric detectors 
Virgo, LIGOs and GEO; however, these progresses should reduce significantly the false 
alarm rate of this search. Further studies with newer and longer data sets are necessary 
to assess quantitatively this issue. Nevertheless, this hybrid network search could be of 
interest for the near future, when only one interferometer will be taking data, a likely 
condition for a large fraction of the time due to the planned instrumental upgrades 
towards enhanced detectors. 
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